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Abstract 

We study J/^p production at fixed-target energies, including the 
colour-octet production mechanisms predicted by NRQCD. In an ear- 
lier paper, we found that the octet components were crucial in under- 
standing pr-integrated forward (xp > 0) cross-sections. Here we make 
a detailed comparison of the theoretical predictions with measured Xp 
distributions from fixed-target experiments. We find that the model 
predictions agree well with data. Taking into account higher orders 
in NRQCD, we argue that cr{xi)/(^{X2) = 0.6, in excellent agree- 
ment with data. Similiar arguments also predict that cj(x)/o'( J/V') is 
roughly equal in photo- and hadro-production. 



1 Introduction 



A large volume of data now exists on the liadro-production of J/ip at fixed 
target energies — both for the total cross section and the Xp distribution. 
In particular, several new experiments have taken high statistics data with 
pion or proton beams. Since a promising new model for heavy- quarkonium 
production is now available, it is useful to perform an analysis in the context 
of this model. This is the purpose of this paper. 

The physics of heavy quarks can be expected to be dominated by its 
(large) mass, m, and hence computable in perturbative QCD. However, for 
quarkonium systems the scales mv (relative momentum) and mv"^ (energy) 
are also important, since the quarks are expected to be non-relativistic in 
the rest frame of the quarkonium. An organised way to include the effects of 
the dimensionless parameter v (relative velocity) is provided by an effective 
low-energy field theory called non-relativistic QCD (NRQCD) [|l|. 

A cutoff of the order of m is introduced into the QCD action to re- 
move all states of momenta larger than m. Their effects are taken into 
account through new couplings, which are local since the excluded states are 
relativistic. Beyond the leading order in 1/m the effective theory is non- 
renormalisable. This cut-off theory is block-diagonalised to decouple the 
heavy-quark degrees of freedom, leaving a non-relativistic Schrodinger field 
theory. Consequently, states of a quark anti-quark pair can be specified by 
the spin, orbital and total angular momentum i^^^^Lj in spectroscopic no- 
tation), a radial quantum number and the colour representation. Gluons 
emitted by quarks can be specified in a multipole expansion. The formalism 
of NRQCD has been successfully applied to first principles lattice computa- 
tions of quarkonium spectra 0. 

In this framework, a cc pair is created by a short- distance process, and it 
binds into a quarkonium state over scales that are longer by powers of 1/v. 
Assuming factorisation, the cross-section for the production of a meson H 
can be written as 



where F„'s are short distance matrix elements and On are local 4-fermion 
operators, of naive dimension describing the long-distance physics (the 
angular brackets denote a vacuum expectation value). The cutoff-dependence 
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of Fn is compensated by that of the long-distance matrix elements. The 
colour index a can run over singlet as well as octet representations. 

The short-distance coefficients F„ have the usual perturbation expansion 
in powers of ag- The scaling of the matrix element {O^i^^^^Lj)) with v is 
determined by the (rotational) tensor representation of the gluonic transition 
connecting the state of the produced pair to the target hadron. Consequently, 
the sum in eq. (|l]) is a double power series expansion in v and Ug- For 
bottomonium states, is small compared to as(m^) and hence colour octet 
contributions are not very significant except when the singlet terms vanish. 
As a result, the expansion is close to the normal perturbative expansion. For 
charmonium states, a numerical coincidence, v"^ ~ as(m^), makes the double 
expansion more complicated. 

The importance of colour octet pieces were first seen in the phe- 
nomenology of P-state charmonium production at large in Tevatron data 
[^]. These processes do not have a consistent description in terms of colour 
singlet operators only However, even for the direct production of S'-states 
such as the J/ip or ip' , where the colour singlet components give the leading 
contribution in f , the inclusion of sub- leading octet states was seen to be 
necessary for phenomenological reasons . This can be understood in terms 
of the numerical coincidence of v and as mentioned earlier. 

Since the long-distance matrix elements are not calculable, and the octet 
matrix elements are not easily written in terms of decays, the Tevatron data 
must be used to fix these unknown parameters As a consequence, quan- 
titative test of this formalism can come only from other experiments. Fixed- 
target experiments provide such tests, because colour-singlet contributions 
are expected to be negligible for photo-production and hadro-production of 
J /if) in appropriate kinematic regions. A new linear combination of octet 
matrix elements appears here. In Ref. ^ this combination has been fixed 
through an analysis of elastic photo-production data. 

In an earlier paper |]TU|, we had considered data on hadro-production of 



J / ijj in proton-nucleon and pion-nucleon collisions at centre of mass energies, 
■\fS, upto about 60 GeV. We found that the data on forward {xp > 0) inte- 
grated cross sections are well described by the inclusion of octet components. 
These are parameter-free predictions, because the values of all the required 
non-perturbative matrix elements are derived from other experiments. In 
this paper, we carry out a more detailed investigation of the Xp and y dis- 
tributions of J/ip production for the same fixed-target energies. These are 
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given in Sections 2 and 3. 

The octet model predictions for J/ip production in fixed-target experi- 



it is claimed 



ments have also been considered in Ref. |TT|, |T2|. In Ref. 
that the ratio of directly produced J/ip to x states is not in agreement with 
data. This is disputed in Ref. |Tl|. We address this question in Section 4. 
Our results are in agreement with those of Ref. [11 



2 Cross Sections 

The Xp distribution for inclusive production of J/ ip can be written as a sum 
over all possible channels in the form 

^= E ^{O^C'^'Lj))BR{H^jm. (2) 

(^^F n={a,S,L,J,H} "'^^ 

A channel is specified by the rotational and colour quantum numbers of the 
cc state, and the colour singlet meson it goes into. For both colour singlet 
and octet channels, the strong interaction effects required to bind a pair into 
a physical meson are specified by the matrix element (O). Subsequently, this 
meson decays into the J/ip with a branching ratio specified by BR. The 
formulae for the channel cross sections |T^, matrix elements and branching 



ratios are collected in Table |l| for all the relevant channels. 

The matrix element appearing in the octet ^Si channel with the gluon 
initiated part of the reaction is the combination 



e = i^or^'s,)) + ^ {ori'Po)) + ^ {or{'P2)) . (3) 

Precisely this combination appears in the cross section for photo-production 
of the J/ip, and has been fitted to the relevant data The two matrix ele- 
ments in the singlet ^Pj state are related to the derivative of the wavefunction 
of xj at the origin by 

orcpj)) = ^-^±^\R'm'. (4) 



It has been fixed by decay rates |14|. Each of the three octet Pj matrix 



elements have been extracted from hadro-production rates at the Tevatron. 
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H 



a 



25+1^ 



7 



da 



BR (%) 



J/ij 
Xo 
Xo 
Xi 

X2 
X2 



L C 3 p 3 p 



2 3 

"^^^ 

2 3 



0.020 ± 0.001 GeV^ 
0.0066 ± 0.0021 GeV^ 
0.32 ± 0.04 GeV^ 

0.0098 ± 0.0013 GeV^ 
0.32 ± 0.04 GeV^ 



100 
100 
0.66 ±0.18 
0.66 ±0.18 
27.3 ± 1.6 
13.5 ± 1.1 
13.5 ± 1.1 



Table 1: The cross sections, non-perturbative matrix elements and branching 
ratios into the J/V^ for each channel (specified by the colour representation 
a, rotational quantum numbers and the primary meson if) relevant to the 
inclusive J/-?/; cross section of eq. (j^). The gluon and quark luminosities are 
defined in eq. @. The value of the matrix element quoted in the first line is 
for the linear combination of eq. (RT). 
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We use the tabulation of [|T5[, along with the relation 



ores,)) = {2j+i){on's^)). (5) 

The value of this matrix element for J = 1 is shown in Table |l|. 

The gluon and quark luminosities, Cg and Cq respectively, are defined in 
terms of the relevant parton distributions in the projectile (P) and target 
(T), through the formula — 

= — =x+gp{x+)x_grr{x^ 



xl + 4r 



xl + 4r 



(6) 

E/ + (P ^ T) , 



where the sum in the definition of Cg is over flavours. We have used the 
notation 



x± = -{^xI + At±Xp). (7) 

We also discuss briefly existing data on rapidity (y) distributions of the J/ip. 
These can be easily obtained from the Xp distribution by multiplying the 
latter with the Jacobian factor \/x'i + 4r for the transformation between x. 



and y. It is clear that the dependence of the cross sections on Xp or y comes 
entirely from these luminosities. 

Before discussing our computations and the comparison with data, we 
would like to note that the colour octet model contains many parameters 
in the form of non-perturbative matrix elements. Our approach here is to 
fix these values from different experiments and use the fixed target hadro- 
production data as a test of the phenomenological viability of the model. 



3 Results 

Before embarking on the computation of cross sections there are the usual 
choices for QCD inputs to be fixed up. From open charm production, the 
limits on the charm quark mass are 1.2 GeV < m < 1.8 GeV 0]. In an 
earlier calculation of the total forward cross section, a{xp > 0) |T^, we had 
varied m between 1.6 GeV and 1.7 GeV to obtain an agreement between the 
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Figure 1: The colour-octet model predictions for integrated forward J/%Ij 
hadroproduction cross sections as a function of the CM energy. For pp colli- 
sions, the two curves for each choice of m and scale Q are for the structure 
functions MRS D— ' and GRV LO. For 7rp coUisions the three sets of structure 
functions are MRS D— ' for proton and SMRS 1 for tt, MRS D— ' for proton 
and SMRS 3 for tt, GRV LO for both. Note that the colour-singlet model 
predictions lie far below the data. 
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data and model computations. The QCD scale, Q, was varied between M 
and 2M (where M is the mass of the J/ip)- For these calculations we had 
used the MRS D— ' and the GRV LO sets of parton densities for the proton, 
and the SMRS 1, SMRS 3 and GRV densities for the pion. All the sets were 



taken from the PDFLIB package JT^- We continue with these choices here. 



The total forward cross sections are shown in Figure |l]. We would like to 
point out the following features — 

• The colour singlet model predictions lie far below the data. The addi- 
tion of the octet channels is necessary for the phenomenology. 

• For pp collisions, all the data lie on the curves shown, with the exception 



of the data from 181. The total forward cross section with a 300 GeV 



proton beam, reported in 1113, is 143 ±17 nb. This is much higher than 



the corresponding number, 71.8 ±9.3, reported by |]T9[ from a 315 GeV 
proton beam. It would be useful to have a third experiment to sort out 
this mismatch. 

• For vrp collisions the data are very scattered. More data would be 
very welcome. The structure functions for vr are also rather badly 
known, and this shows up as a bigger spread in the model prediction. 
Further experiments on high-mass lepton pair production with pion 
beams could reduce this uncertainty. 

• The data at the highest energies (ISR) seem to prefer a smaller value 
of m. We discuss this point later. 

In Figure ^, we compare the model with data on Xp distributions from 



three different pp experiments |T8|, EDI |21|. There is extremely good agree- 



ment for proton beam energy of 800 GeV over the 4 orders of magnitude 
spanned by the data pi|. Since the data sample contains 2.45 x 10^ events. 



the quoted errors are very small and the agreement is very significant. Data 
at other beam energies are also well described. The Xp distribution at 300 
GeV energy [|1^] has been normalised by the ratio of their measured value 
of a{xp > 0) and the model prediction. There is no statistically significant 
evidence of a systematic deviation from the model predictions at large Xp. 
At large values of Xp the 300 GeV data may lie slightly below the model 
(but with statistically insignificant deviation). On the other hand, data in 
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the last Xp bin for the 800 GeV set hes a httle above the prediction. Note 
also that only the statistical errors on the data are shown in the figures. 

Figure ^ shows the Xp distributions in n—p collisions. The data are taken 
from three different experiments [0, |2^, At the highest beam energy 



of 515 GeV, the data agrees extremely well with the model predictions both 
in magnitude and shape. At lower beam energies the shape is reproduced 
extremely well, although the agreement in magnitude is less perfect. It is not 
possible to decide whether this is due to nuclear effects or the pion structure 
functions, or even experimental uncertainties in the normalisation (see Figure 
|1|). More data would help in deciding between these possibilities. 

There is a small body of literature on cross sections and distributions in 
pp collisions. The data on total forward cross sections is even more scattered 
than the pion data. Also, the number of J/ ip events are rather low, and hence 
lead to larger errors. For these reasons we have not attempted a quantitative 
study of this data. 

Since a large amount of data exists for the differential cross section da/dy 
at y = in pp collisions, we have also tested the model against this. The 
data goes from values of near threshold upto ISR energies. The compar- 
ison of model and data are shown in Figure ^ It is clear that the data on 
the rapidity distribution is more scattered than that on total forward cross 
sections. In several cases experiments at the same value of -sfS yield values 
for the distribution which are not in agreement even when errors are taken 
into account. Nevertheless, a reasonable overall agreement of the model with 
the data can be seen. 

It is interesting to note, however, that the low data prefer a larger 
value of the charm quark mass, m, and a smaller scale Q than the ISR 
data. As shown in the figure, this trend holds for both the total forward 
cross section and the rapidity distribution. A similiar observation used to be 
made for single inclusive charm production, before higher order perturbative 
calculations resolved this problem The essential physics is that higher 
orders are important near threshold because an off-shell heavy quark can be 
thrown on-shell by soft radiation. It seems reasonable to expect that similiar 
effects are at work here. Further study would call for inclusion of higher 
order perturbative effects. 

Very naively one may have thought of tuning m and Q in order to produce 
the best results. However, since such a tuning would depend on v^, it is 
clearly not the correct way to proceed. In fact, as is now well-known, such 
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Figure 2: The colour-octet model predictions for the cross section differential 
in Xp for J/ip production in pp collisions for different beam energies. The two 
curves are for the parton density sets MRS D— ' and GRV LO and m = 1.7 
GeV and Q^M. 
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Figure 3: The colour-octet model predictions for the cross section differential 
in Xp for J/ ip production in np collisions for different beam energies. The two 
curves are for the parton density sets MRS D— ' and GRV LO and m = 1.7 
GeV and Q = M. For visibility, the 7r+ and vr^ data for beam energy 300 
GeV have been offset equally to either side of the true Xp value. This is also 
done for the Cu and Be data at 125 GeV. 
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Figure 4: a{xp > 0) for J/ip production in pp collisions (upper panel) and the 
differential cross section Bda/dy at y=0 (lower panel), shown as a function 
of \/S. For the model calculation we have assumed that the branching ratio 
for the decay of the J/ijj into leptons (B) is 5.97%. 
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Source of uncertainty 


pp 


np 


Parton density sets 


20% 


29% 


Scale uncertainties 


91% 


96% 


Non-perturbative matrix elements 


5% 


5% 



Table 2: Estimates of the magnitude of various theoretical uncertainties in 
the colour octet model for J/^p production at = 19.4 GeV. These are 
obtained from a calculation of da/dxp at Xp = 0. 



sensitive dependence on the scale is an artifact of leading order perturbative 
QCD calculations. Once higher orders are computed, it is possible to remove 
this strong dependence on Q and then fix the best values for m |jT6|. Hence, 
we may take the joint dependence of the cross sections on m and Q as an 
estimate of the importance of higher order perturbative QCD effects. 

The magnitudes of the various theoretical uncertainties are summarised in 
Table ^ The errors on the individual non-perturbative matrix elements are 
quoted in Table |I|. The dominant channel in eq. (Q) is the gluon fusion part of 
the octet J/ ip (the other channels are suppressed either by a branching ratio 
or by Cg/Cg). As a result, the corresponding matrix element dominates this 
source of uncertainty. The error quoted in the table above is a 1-a limit. Even 
if it is replaced by a 3-a limit, the uncertainty due to the matrix elements 
remains the smallest source of error. Scale uncertainties in the cross section 
are a factor of two. As discussed in the previous paragraph, these may be 
taken as a rough estimate of the importance of higher order perturbative 
QCD effects. 

Apart from the model discussed in this paper, another, called the semi- 
local duality model (SLD) [|^, is also used to describe charmonium cross 
sections. This model fixes cross sections for individual charmonium states 
by comparison with data on the total cross section at some value of v^- 
Once this is done, the results for dependence and similiar to 

those reported here. A distinction can perhaps be made with accurate data 
at large Xp, since SLD predicts a somewhat harder Xp distribution. Other 
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Figure 5: i?, the fraction of da/dxp as a function of Xp in to the channels 
(a) gg J/iIj, (b) qq Jji^), (c) Xo, (d) Xi and (e) X2 production. 

distinctions are possible [|11[], and will be discussed elsewhere. 

4 Discussion 

Next we turn to the ratio o"(x)/o"(^/^)- All the hard cross sections in Table 
m are taken to order as- The non-perturbative matrix elements in the Xj 
channels are all of leading order in the velocity — f ^, and there are no con- 
tributions at order w^. The octet terms for Xj are suppressed relative to the 
singlet terms by Cqj Cg. The contribution to direct J/ip production starts at 
order f^. Although sub-leading in f , it happens to be the channel with the 
largest contribution (see Figure ^ since it is gluon initiated and hence comes 
with a factor of Cg. The xi state, which has the largest branching fraction 
into J/ifj has no singlet contribution and hence is also by Cq/ Cg. 

From this discussion it is clear that the relative magnitudes of the parton 
luminosities, Cg and Cq, are important in determining the dominant produc- 
tion channels, and may offset the power counting in v to some extent. At 
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order there are no new contributions to direct J/ip production, but xi 
is produced through gluon initial states. Therefore, we have to consider all 
production channels upto order to get the x to direct J/ip ratio. Beyond 
the order terms, power counting in v can be used to truncate the NRQCD 
series at any desired level of accuracy, since all the dominant Cg type terms 
are accounted for at order v^. 

Enumerating all possible processes at this order, we find that a certain 
linear combination, of (Of'^So)), (Of (=^Po,2)) and (Of (^DjO) is re- 
quired for the hadro-production cross section ||2^. These arise equally from 
quark and gluon initiated processes and hence the same linear combination 
$ is associated with both Cg and Cq. Analysis of the diagrams shows that 
precisely the same linear combination appears also in the photo-production 
cross section. 

Two interesting consequences now emerge as predictions of the NRQCD 
approach to quarkonium production — 

1. There are three families of matrix elements involved in Each scales 
as (2 J -|- 1) (with order f ^ corrections) [0] where J is the spin of the x 
state. This gives — 

^ ^ I I (NRQCD), 
^^^^^ [ 0.62 ±0.18 ±0.09 (experiments). 

Within NRQCD, the predicted ratio does not depend on whether the 
mesons are photo- or hadro-produced. Nor does it depend on v^- 
The measurement quoted above comes from a global analysis of hadro- 
production experiments [^. There are no studies of this ratio in photo- 
production. We would like to urge experimentalists to perform such a 
measurement. 

2. In photo-production the ratio of x and direct J/ip production is deter- 
mined by the the ratio of the terms in $ and G (eq. |^). As a result 
this ratio is independent of y/S. In hadro-production this ratio is al- 
tered a little by the contribution of the colour- singlet production of Xo 
and X2- The ratio increases by about 0.05 (see Figure ||). A small 

dependence is also to be expected in hadro-production due to the 
\^ dependence of Cq/Cg. This may cause the ratio to be slightly 
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larger at small values of v^- A detailed numerical treatement will 



be presented elsewhere [27|. Data on this ratio from hadro-production 



yields the value 0.43 ± 0.04 ± 0.04 pSI, almost independent of VS. In 



photo-production there is only the limit (xi + X2)/{J/'^) < 34% (90% 
CL) [^. Data on this ratio from photo-production experiments would 
constitute an important test of the NRQCD based colour-octet model. 

The introduction of these order terms does not change the inclusive 
J/ip cross section because the value of G extracted from photo-production 
has to decrease due to the inclusion of $. The shape of the Xp distribution 
comes from the gluon luminosity. This is clear from Figure ^, since the ratio 
of the direct to total J/ if) differential cross section is constant. Since the order 
f ^ contribution also comes from a gluon initiated process, and carries a factor 
Cg, the shape of the Xp distribution remains unchanged. The net result of 
adding all contributions to the cross section upto order is to change only 
the X to direct J/ip ratio. A detailed study will be published later p7[ . 



In summary, the colour-octet model based on NRQCD and an assumption 
about factorisation, written to order agV^ gives a good phenomenological de- 
scription of the total forward J / ip cross section and the Xp and y distributions 
in hadroproduction at fixed target energies. A small discrepancy is found in 
extrapolating the fixed target results to ISR energies, but these are of a kind 
that have been earlier seen in single inclusive charm production and removed 
there by going to higher orders in as- There is a discrepancy in the ratio of 
X to direct J/ip production at order v^. We have argued that this problem 
can be solved without changing the results on the cross section by extend- 
ing the calculations to order asV^. Inclusion of these order terms leads 
to two predictions — that for the Xi/X2 ratio is borne out by data, and the 
other, the rough equality of the xl J 1'^ ratio in photo- and hadro-production 
can easily be tested in experiments. We beheve that NRQCD provides a 
phenomenologically viable model for quarkonium production. Control of the 
next orders in the expansions in v and as are, however, necessary. This work 



is now in progress p7 |. 
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